***Datafolha 2005
use "C:\Users\mauricioizumi\Desktop\notebook\BPSR\2020-0009\02528.dta", clear

***ESCOLARIDADE
tab escola [iweight=pesoe]
gen escolaridade = escola
replace escolaridade = 5 if escola == 6
replace escolaridade = 7 if escola == 8
replace escolaridade = 3 if escola == 4
tab escolaridade [iweight=pesoe]


***RENDA
tab rendaf [iweight=pesoe]
gen renda = rendaf
replace renda =4 if rendaf == 5 | rendaf == 6 | rendaf == 7
replace renda=. if rendaf == 99
tab renda [iweight=pesoe]



***Tabela 1
tab p14 [iweight=pesoe] 
gen avaliacao_lula_categorica = p14 
replace avaliacao_lula_categorica =. if p14 == 99
replace avaliacao_lula_categorica =1 if p14 == 2
replace avaliacao_lula_categorica =4 if p14 == 5
tab avaliacao_lula_categorica [iweight=pesoe]


tab PARTIDO [iweight=pesoe]
gen IP = PARTIDO
replace IP = 10 if PARTIDO == 1 | PARTIDO == 3 | PARTIDO == 4 | PARTIDO == 5 | PARTIDO == 7 | PARTIDO == 8 | PARTIDO == 9 | PARTIDO == 10 | PARTIDO == 11
replace IP =0 if PARTIDO == 12 | PARTIDO == 13
tab IP [iweight=pesoe]


gen voto_1t=.
replace voto_1t=1 if p11a == 4
replace voto_1t=0 if p11a == 7
replace voto_1t=2 if p11a == 1 | p11a == 2 | p11a == 3 | p11a == 6
tab voto_1t [iweight=pesoe]




****Tabela 2 
**Recebe PBF - entrevistado
tab p28aa [iweight=pesoe]

***Somando Entrevistado e Fam�lia - PBF
gen recebe_pbf = p28aa
replace recebe_pbf = 1 if p28aa == 2
replace recebe_pbf = 1 if p28ab == 2 
replace recebe_pbf = 0 if p28aa == 3
replace recebe_pbf = 0 if p28aa == 96
replace recebe_pbf = 0 if p28ab == 3
replace recebe_pbf = 0 if p28ac == 3
tab recebe_pbf [iweight=pesoe]


**Recebe Bolsa Escola - entrevistado
tab p28ba [iweight=pesoe]


***Somando Entrevistado e Fam�lia - Bolsa Escola
gen recebe_escola = p28ba
replace recebe_escola = 1 if p28ba == 2
replace recebe_escola = 1 if p28bb == 2 
replace recebe_escola = 0 if p28ba == 3
replace recebe_escola = 0 if p28ba == 96
replace recebe_escola = 0 if p28bb == 3
replace recebe_escola = 0 if p28bc == 3
tab recebe_escola [iweight=pesoe]




*****VD: Avalia��o
mlogit avaliacao_lula_categorica i.recebe_pbf i.recebe_escola ib4.regiao i.renda i.SEXO i.escolaridade i.IDADE  i.metrop [pweight=pesoe], baseoutcome (4) rrr
estimates store Aprova
**Tabela 3
mchange recebe_pbf recebe_escola, stats(change from to pvalue ci) brief 


****VD: IP
mlogit IP i.recebe_pbf i.recebe_escola ib4.regiao i.renda i.SEXO i.escolaridade i.IDADE  i.metrop [pweight=pesoe], baseoutcome (0) rrr
estimates store IP
**Tabela 3
mchange recebe_pbf recebe_escola, stat(change from to pvalue ci) brief 



****VD: Voto primeiro turno
mlogit voto_1t i.recebe_pbf i.recebe_escola ib4.regiao i.renda i.SEXO i.escolaridade i.IDADE  i.metrop [pweight=pesoe], baseoutcome (0) rrr
estimates store Voto
**Tabela 3
mchange recebe_pbf recebe_escola, stat(change from to pvalue ci) brief 



**Gr�fico - 1
coefplot (Aprova, label (VD: Aprova Governo Lula (em relação a quem não aprova) N=2444)) (IP, label (VD: Identificação com PT (em relação a quem não tem identificação) N=2476)) (Voto, label (VD: Intenção de Voto em Lula (em relação a Alckmin) N=1854)), keep(*recebe_pbf *recebe_escola) eform xline(1) legend(rows(1))  xtitle(Razões de Risco Relativo) title (Impacto de ser Beneficiário sobre Opiniões Políticas em 2005 - 95% IC - Multinomial logística)




***ESEB 2018
use "C:\Users\mauricioizumi\Desktop\notebook\BPSR\2020-0009\04622.dta", clear


****Recebe PBF
tab P29
gen recebe_pbf =.
replace recebe_pbf =1 if P29==1
replace recebe_pbf =0 if P29==2
tab recebe_pbf

tab D1A_FAIXAID
gen idade=.
replace idade=1 if D1A_FAIXAID == 1 | D1A_FAIXAID == 2
replace idade=2 if D1A_FAIXAID == 3
replace idade=3 if D1A_FAIXAID == 4
replace idade=4 if D1A_FAIXAID == 5
replace idade=5 if D1A_FAIXAID == 6 | D1A_FAIXAID == 7
tab idade 


***G�nero
tab  D2_SEXO
gen sexo = D2_SEXO


***Escolaridade
tab D3_ESCOLA 
gen escola=.
replace escola=1 if D3_ESCOLA == 0 | D3_ESCOLA == 1
replace escola=2 if D3_ESCOLA == 2 | D3_ESCOLA == 3
replace escola=3 if D3_ESCOLA == 4 | D3_ESCOLA == 5
replace escola=4 if D3_ESCOLA == 6 | D3_ESCOLA == 7
replace escola=5 if D3_ESCOLA == 8 | D3_ESCOLA == 9
tab escola


***Renda
tab D9
gen renda_parcial=. 
replace renda_parcial=1 if D9 <= 954 /* At� 1 SM */
replace renda_parcial=2 if D9 > 954 & D9 <= 1908 /* De 1 a 2 SM */
replace renda_parcial=3 if D9 > 1908 & D9 <= 4770 /* De 2 a 5 SM */
replace renda_parcial=4 if D9 > 4770 /* Acima de 5 SM */
replace renda_parcial=. if D9 == 9999998 | D9 == 9999999

gen renda=renda_parcial
replace renda=1 if D9A==1
replace renda=2 if D9A==2
replace renda=3 if D9A==3
replace renda=4 if D9A==4 | D9A==5 | D9A==6 | D9A==7
tab renda





**Tabela 4
gen avaliacao_temer=.
replace avaliacao_temer=0 if Q9== 1 | Q9== 2
replace avaliacao_temer=1 if Q9== 3
replace avaliacao_temer=2 if Q9== 5 | Q9== 5 
tab avaliacao_temer

tab Q10A
tab Q10B
gen IP_PT =2
replace IP_PT =1 if Q10B == 13
replace IP_PT =0 if Q10A == 5 
replace IP_PT =. if Q10B == 97 | Q10B == 98 | Q10A == 7 | Q10A == 8
tab IP_PT


tab P1
tab P1A
gen Anti_PT =2
replace Anti_PT =1 if P1A == 13
replace Anti_PT =0 if P1 == 2 
replace Anti_PT =. if P1A == 97 | P1A == 98 | P1A == 99 | P1 == 8 | P1 == 9 
tab Anti_PT

gen voto_1t=.
replace voto_1t=0 if Q12P1_B == 5
replace voto_1t=1 if Q12P1_B == 9
replace voto_1t=2 if Q12P1_B == 1 | Q12P1_B == 2 | Q12P1_B == 3 | Q12P1_B == 4 | Q12P1_B == 6 | Q12P1_B == 7 | Q12P1_B == 8 | Q12P1_B == 10 | Q12P1_B == 11 | Q12P1_B == 12 | Q12P1_B == 13 | Q12P1_B == 14
tab voto_1t


****VD: avalia��o Temer
mlogit avaliacao_temer i.recebe_pbf i.renda i.escola i.sexo i.REG i.idade, cluster(UF) rrr baseoutcome(2)
estimates store Aprova
**Tabela 5
mchange recebe_pbf, stat(change from to pvalue ci) brief

****VD: IP
mlogit IP_PT i.recebe_pbf i.renda i.escola i.sexo i.REG i.idade, cluster(UF) rrr 
estimates store IP
**Tabela 5
mchange recebe_pbf, stat(change from to pvalue ci) brief

***VD: anti-PT
mlogit Anti_PT i.recebe_pbf i.renda i.escola i.sexo i.REG i.idade, cluster(UF) rrr baseoutcome(0) /* Sem antipatia � baseline */
estimates store anti_IP
**Tabela 5
mchange recebe_pbf, stat(change from to pvalue ci) brief

****VD: voto primeiro turno
mlogit voto_1t i.recebe_pbf i.renda i.escola i.sexo i.REG i.idade, rrr cluster(UF) baseoutcome(1)
estimates store Voto
**Tabela 5
mchange recebe_pbf, stat(change from to pvalue ci) brief

**Gr�fico - 2
coefplot (Aprova, label (VD: Aprova Governo Temer (em rela��o a quem n�o aprova) N=1750)) (IP, label (VD: Identifica��o com PT (em rela��o a quem n�o tem identifica��o) N=2076)) (anti_IP, label (VD: Rejeita PT (em rela��o a quem n�o rejeita nenhum partido) =2100)) (Voto, label (VD: Voto em Haddad (em rela��o a Bolsonaro) N=1610)), keep(*recebe_pbf) eform xline(1) legend(rows(1))  xtitle(Raz�es Risco Relativo) title (Impacto de ser Benefici�rio sobre Opini�es Pol�ticas em 2018 - 95% IC -  Multinomial log�stica)


